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We solve the coupled gap equations for the a- and 7r-bands of MgB2 in the vortex state and cal- 
culate the resulting field dependencies of the specific heat coefficient 7 and the thermal conductivity 
k . The crucial parameters of the theory are the interband pairing interaction A T(T and the ratio 
s — ^a/^TT of the coherence lengths. For reasonably small A WCT and s, the small gap decreases 
with increasing magnetic field H much faster than the large gap A CT . This gives rise to the observed 
rapid increase of 7^ and k-k for small fields while j a and n a exhibit conventional field dependencies. 
Inclusion of intraband impurity scattering yields fairly good agreement with experiments for applied 
fields along the c axis. 



Evidence for the for the existence of two superconducting gaps in MgB2 1 is provided by the rapid rise of the specific 
heat coefficient j s {H) 2 and the thermal conductivity k s (H) 3 at very low fields. These measured field dependencies 
can be explained qualitatively by assuming two independent bands where the large s-wave pairing gap A CT = Ai 
is associated with the two-dimensional cr-band and the small s-wave gap A T = A2 is associated with the three- 
dimensional 7r-band. 4 The steep rise of js(H) and k s (H) can be explained qualitatively by assuming that the "virtual" 
upper critical field for A v is much smaller than that of A^ ? In the present paper we improve the theory of Ref. 4 
by taking into account the interband pairing interaction while neglecting the interband impurity scattering 5 which 
has been shown to be small. 6 We first have to solve the two-band gap equations in the presence of the vortex lattices 
produced by a magnetic field. Generalization of the linearized gap equations near the upper critical field H C 2 7 to all 
averaged fields H between H C 2 and yields, instead of the single gap equation of Ref. 4, the following coupled gap 
equations for the gaps Ai at T = : 



dwB:i(u,Aj,A/vj), (i = l,2) 



(1) 



Here, the Ay are the intra- and interband pairing interactions multiplied by the densities of states iVj(O) , and the Bj 
are the spectral functions of the anomalous propagators for the Abrikosov vortex lattice: 4 



B, = Rc 



-i ypK 2 A ri (A/ Vi) w(zi) 



{1 + SA^A/^Hl + i^ Zl w( Zl )]} 1/2 



(2) 



Zi = 2(oj + A/vi ; A= (2e£T) _1/a ; ll = T t A l 



(3) 



A ri = Ai/A ; A = 1 - 2 (7iA/t>i) w{ Zi ) . 



(4) 



The Ti are the normal state impurity scattering rates, the Ai are the field-dependent densities of states, the Vi are 
the Fermi velocities perperdicular to the field, and the function w(z) is defined in Ref. 4. For band 2 (the 7r-band) we 
assume a spherical Fermi surface. Then vi is replaced by vi sin 9 and B2 and A are averaged over the polar angle 9 
with respect to the direction of H . For brevity we omit here, and in the following, the terms containing s'm9 and the 
integrations over d9 sin 9 from to 7r/2. For band 1 (the cr-band) and the field along the c axis, 9 — tt/2 and thus 
sin 9=1. 

In the limit H —> the vortex lattice constant a = (2/7r) 1 / 2 A tends to infinity. Making use of the asymptotic 
expansion w(z) ~ if ' ypn z , the gap equations, Eq.(l), become 



A l0 = J2 Xi - 



duj Re 



(- 2 ~ A| c 



where the A^o are the gap values at T = in zero field. Here we have made use of the relation 

(u> + iji ) I A° ri = (uj + i-f i0 ) Ao/Aio = uj/A i0 



(5) 



(6) 
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which can be derived with the help of the expression for Di in Eq.(4). However, this relation only holds in the absence 
of interband impurity scattering. 5 

With the help of the Abrikosov parameters, denoted by ft , we now express A/vt in terms of the reduced field 
h = H/H C 2 and the zero-field gap A i0 : 4 

A/vi = (6ft h A 2 )~ 1/2 ; A/v 2 - a (6ft h A^)" 1 / 2 ; s = ( Ul /« 2 )(A 20 /A 10 ) = W&o • (7) 

Employing these relations we can express the gaps Aj and the scattering rates Ti in Eq.(2) by their ratios with respect 
to the Aio , and we can convert the integrations over u> in Eq.(l) to integrations over the new variables f2 = w/A^o • 
We then divide Eq.(l) by Aj and Eq.(5) by A^o and subtract the latter from the former. In this way we obtain two 
coupled equations for the two unknown functions x\{h) and x 2 (h) for given values of the parameters Ay , Si, ft , r , 
and s . The quantities , Si , and r are defined by 

x, = A,/A 40 ; Si = r,/A 40 ; r = A 20 /A 10 . (8) 

In Fig. (1) we show the reduced gap functions x\(h) and x 2 (h) for 3 sets of parameter values: I) An = 1 , A 22 = 0.28 , 
A21 = 0.17, A12 = 0.23, ft = 1.15, ft = 1.58, s = 1/6; la) A y - as in set I, s = 1/4, ft = 1.21, ft = 1.58; II) 
An = 1 1 A 22 = 0.45 , A 2 i = 0.16 , Ai 2 = 0.21 , ft = 1.16 , ft = 1.58 , s = 1/6 . The gap ratio has been taken to be 
r = 1/3 and the reduced impurity scattering rates Si — I^/A^o are S\ — 0.5 and S 2 = 0.8 . The A-matrix elements in 
I) and II) have been obtained from band structure calculations (Refs. 8 and 9, see the review, Ref. 10). The ratio of 
Fermi velocities is about vi/v 2 = 0.54, 11 and the gap ratio r ranges between about 1/3 and 0.44 12 ' 13 which yields a 
range of the ratio s of coherence lengths (see Eq.(7)) between 0.18 and 0.24. While the function Xi(h) ~ (1 — h) 1 ^ 2 
is rather independent of the choice of parameters, the function x 2 (h) depends sensitively on the values of Ajj and s . 
For vanishing interband coupling A 2 i we obtain approximately x 2 (h) ~ (1 — h/ s 2 ) 1 / 2 which goes to zero at a smaller 
effective upper critical field H£ 2 — s 2 H c2 where s ~ ^10/^20 (see Eq.(7)). Thus we see that H£ 2 corresponds to the 
"virtual" upper critical field for the 7r-band which was introduced in Ref. 2 as the field above which the overlap of the 
vortex cores with large radius £ 2 o (see the STS measurements of Ref. 12 ) drives the majority of the 7r-band electrons 
normal. 

In Fig. (2) we have plotted our results for the zero energy density of states Ai(u> = 0) obtained from the expression 
for Ai (given by Eq.(2) with the numerator set equal to 1) by inserting the previously calculated gap ratios Xj(/i) 
together with Eqs.(7) for the functions A(h)/vi . A 2 is obtained by averaging A 2 {9) over the polar angle 9 . We note 
that it is important to calculate the impurity scattering rates 7i = TiAi(h) self-consistently. One sees from Fig. (2) 
that A 2 (h) rises steeply for small fields h and then becomes almost constant above h <~ 0.2 . The slope at h = and 
the downward curvature for low fields increase as s is decreased from 1/4 to 1/6. The function A\{h) is very similar 
to the function obtained previously for a single band. 4 The initial steep rise of A 2 (h) for s = 1/6 qualitatively fits the 
data points for the contribution of the 7r-band to the specific heat coefficient -f(H) 2 The function A\(h) corresponds 
to the straight line assumed in Ref. 2 for the cr-band contribution to ^(H) for fields applied along the c axis. 

We turn now to the calculation of the in- plane electronic thermal conductivity K s (h) which is given at T = by 
the expression in Ref. 14. Again it is important to take into account the renormalization of the gap by the function 
Di in the presence of impurity scattering (see Eq.(13) of Ref. 4). By inserting the functions Xi(h) obtained from the 
gap equations, Eq.(l), and the functions A(h)/vi from Eq.(7) into the expresions for the ratios n s /K n we obtain, for 
applied fields H along the c axis, the thermal conductivity ratios n si {h)/ n ni shown in Fig. 3. It should be noted that 
the 7r-band contribution K s2 (h)/n n2 has been obtained as an average over the polar angle 9 by including the factor 
(3/2) sin 2 9 which arises from the square of the group velocity in the ab plane. The curve for the c-band conductivity 
n s i(h)/n n i turns out to be very similar to the curve obtained in Ref. 4 for a single band with the same impurity 
scattering rate 5\ — 0.5 . The curve n s2 (h)/n n2 for the 7r-band contribution with S 2 = 0.8 rises almost linearly with 
h where the slope near h = and the downward curvature for low fields increase as s is decreased from 1/4 to 1/6. 
For applied fields perpendicular to the c axis, the measured thermal conductivity K e first rises steeply for small fields 
and then saturates while, for fields along the c axis, it exhibits an upward curvature towards H c2 . 3 These different 
behaviors have been explained by separating the individual contributions of the it- and c-bands. Then rises steeply 
with H and approximately attains its normal-state value at a small field and n a first rises very slowly and then curves 
upward towards H c2 . 3 These experimental curves are similar to our results for s = 1/6 shown in Fig. (3). 

We now briefly discuss the parameter values and approximations that have been used to derive our results. We have 
seen that the rapid increase of the specific heat coefficient j s (h)/j n and the thermal conductivity ratio n s {h)/n n with 
increasing field h = H/H c2 is due mainly to the 7r-band contribution. The reason is that the gap A 2 (/i) associated with 
the 7r-band almost closes at the so-called "virtual" upper critical field 2 H* 2 <~ s 2 H c2 because the ratio s — £ w /£ 20 
of the coherence lengths of the a- and the 7r-bands is much smaller than 1. For the most important parameters 
entering our gap equations, s and r = A 20 /Ai , we have used the values r — 1/3 and s = 1/4, and 1/6 which are 
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based on various experiments on MgB2 . 2 ' 11-13 The field dependence of the gap ratio x\(h) = A1/A10, and thus 
of the contributions to 7si/7ni and K s \/K n i arising from the cx-band, are nearly the same as those obtained for an 
independent single cr-band, 4 indicating that the effect of the interband coupling A12 is rather small. However, the field 
dependence of x 2 (h) = A 2 /A 2 o differs substantially from that for the independent single 7r-band with an effective 
upper critical field h™ 2 = s 2 as can be seen from Fig.(l). This is because X2(h) is non-zero between h* 2 and h = 1 due 
to the effect of the interband coupling A21 . This shows that superconductivity survives even though the vortex cores 
for the 7r-band with giant radius £20 12 start to overlap for h > h^ 2 ■ As can be seen in Fig.(l), the curve for X2(h) 
is very sensitive to the values of Ajj and s . We find that the experimental contributions to 7 and n arising from the 
7r-band 2 ' 3 can be fitted by taking the \{j given by band structure calculations. 8-10 The other crucial parameter values 
needed to obtain good fits of the data are r = 1/3 and s = 1/6 which lie in the ranges obtained from experiment. 
The other parameter values used in our numerical calculations are the reduced impurity scattering rates Si = 0.5 and 
S 2 = 0.8 which have been estimated from the relevant experiments. It turns out that even for these moderately large 
impurity scattering rates it is very important to take into account the renormalization of the gap (see Eq.(4)) which 
leads to a large reduction of the effect of impurity scattering in comparison to that calculated without the function 
D. It is also important that the calculation of the scattering rate 7* = TiAi(h) in the Born limit be carried out 
self-consistently together with the calculation of the zero-energy density of states Ai{h) . This yields Aj(0) = , as 
it should. The shape of the Fermi surface (FS) and the direction of the applied field play an important role because 
the spectral functions Bi in Eq.(2) have to be averaged over the corresponding FS where the velocity Vi denotes the 
component Vi±(p) perpendicular to H . For the spherical FS we have assumed for the 7r-band, v±(p) = vsin0 , where 
9 = l(p,H) . In the limit 9^0 the function B approaches the BCS spectral function of the anomalous propagator. 
We find that the results for the averages over the polar angle 9 do not differ significantly from the results obtained 
by setting 9 = n/2 . This means that the quasiparticles moving perpendicular to the vortex axes yield the dominant 
contributions to 7 and k . We have approximated the 7r-band FS by a sphere whereas the Ay have been calculated 
for the actual FS. This actual FS can be modeled by a half-torus 7 which yields, with the A-matrix of Ref. 8 and small 
s, results 15 which agree qualitatively with ours shown in Figs. 1 and 2 for s = 1/6. Finally it should be pointed 
out that we have employed the Abrikosov ground state of the vortex lattice although, in particular at lower fields, 
a Landau-level expansion or a variational expression 7 is needed to describe the distorted vortex lattice. The results 
for 7 a (ft)/7 n and K s {h)/n n shown in Figs. (2) and (3) for the a- and 7r-bands should still be added by weighting them 
with the corresponding density of states. 

In conclusion we can say that our two-band theory for the vortex state in MgB 2 can satisfactorily account for 
the observed field dependence of the specific heat coefficient 7 and the thermal conductivity k . The small gap A2 
associated with the 7r-band decreases with increasing field H much faster than the large d-band gap Ai which shows 
conventional field dependence. This gives rise to the rapid increase of 7 and k at small fields. Due to a small interband 
pairing interaction A21 , the gap A2 remains finite even in the field region where the large 7r-band vortex cores of radius 
£20 overlap. This leads to smooth evolution of the 7r-band contributions to 7 and k to their normal state values near 
a "virtual" upper critical field H£ 2 — s 2 H c2 which is much smaller than H c2 because the ratio s = ^10/^20 is much 
smaller than 1. The Fermi surface topology and the impurity scattering have relatively small influence on the field 
dependence of 7 and n for fields applied along the c axis. 

We thank T. Dahm for valuable discussions. 
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Fig. 1 




FIG. 1. Reduced gaps xi(h) = A1/A10 for the cr-band and x 2 (h) = A2/A20 for the 7r-band vs reduced magnetic field 
h — H/H c2 for applied fields along the c axis. The upper curves are for x\ and the lower curves for x 2 ■ The pairing interaction 
matrices Xij have the values given in Ref. 8 (our parameter sets I, s = 1/6, and la, s = 1/4) and in Ref. 9 (set II, s — 1/6) 
where s is the ratio of coherence lengths s = ^10/^20 • The curves correspond to the parameter sets la, I, and II, from top to 
bottom. The reduced impurity scattering rates 5t are 5i = 0.5 for the cr-band and 8 2 = 0.8 for the 7r-band. 
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Fig. 2 




H / H C2 



FIG. 2. Specific heat coefficients, or zero-energy densities of states, j s i/~/ n i = M , (i — 1, 2) vs h for the parameter values of 
Fig. (I). The lower curves are for the <7-band, i = 1 , and the upper curves are for the 7r-band, i = 2 , for parameter sets la, I, 
and II, from bottom to top. 
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Fig. 3 




H / H C2 



FIG. 3. Reduced electronic thermal conductivities, K si /K n i , (i = 1,2) vs h for the parameter values of Fig.(l). The lower 
curves are for the cr-band, i = 1 , and the upper curves are for the 7r-band, i — 2 , for parameter sets la, I, and II, from bottom 
to top. 
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